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The electronic environment causes decoherence and dissipation of the collective surface plasmon 
excitation in metallic nanoparticles. We show that the coupling to the electronic environment influ- 
ences the width and the position of the surface plasmon resonance. A redshift with respect to the 
classical Mie frequency appears in addition to the one caused by the spill-out of the electronic density 
outside the nanoparticle. We characterize the spill-out effect by means of a semiclassical expansion 
and obtain its dependence on temperature and the size of the nanoparticle. We demonstrate that 
both, the spill-out and the environment-induced shift are necessary to explain the experimentally 
observed frequencies and confirm our findings by time-dependent local density approximation cal- 
culations of the resonance frequency. The size and temperature dependence of the environmental 
influence results in a qualitative agreement with pump-probe spectroscopic measurements of the 
differential light transmission. 

PACS numbers: 78.67.Bf, 73.20.Mf, 71.45.Gm, 31.15.Gy 



I. INTRODUCTION 

One of the most prominent features of a metal- 
lic nanoparticle subject to an external driving field is 
a collective electronic excitation, the so-called surface 
plasmoni^^ Since the first spectroscopic measurement 
of the related resonance in the absorption cross section 
of free sodium clusters^ much progress has been made 
in the characterization of this collective resonance, both 
experimentally 1 ! 5 ! 6 ! 7 ! 8 : 9 and theoretically. 2 ! 10 ! 11 ! 12 ! 13 ! 14 
The proposed application of metallic nanoparticles^ or 
nanocrystalsi^ as markers in biological systems such as 
cells or neurons renders crucial the understanding of their 
optical properties. 

The first experiments have been made on ensem- 
bles of nanoparticles, where the inhomogeneous broad- 
ening of the resonance resulting from the size depen- 
dence of the resonance frequency masks the homogeneous 
linewidthiiLiSii 9 . In order to gain detailed information 
on the collective resonance, considerable effort has lately 
been devoted to the measurement of single-cluster optical 
propertiesi^iSiiSSiS 3 ^ The possibility of overcoming the 
inhomogeneous broadening resulted in a renewed interest 
in the theory of the optical response of metallic clusters. 

From a fundamental point of view, surface plasmons 
appear as interesting resonances to study given the vari- 
ous languages that we can use for their description, which 
are associated with different physical images. At the clas- 
sical level, a nanoparticle can be considered as a metallic 
sphere of radius a described by a Drude dielectric func- 
tion e(u>) = 1 — u>p/u)(uj + vyi), where u> p — y/ '47m e e 2 /m e 
is the plasma frequency, 7; the relaxation or collision 
time, while e, m e and n c stand for the electron charge, 
mass and bulk density, respectively. Classical electro- 
magnetic theory for a sphere in vacuum yields a reso- 



nance at the Mie frequency um = t^p/y^i 1 ' 2 ' 3 

At the quantum level, linear response theory for an 
electron gas confined by a positive jellium background 
yields a resonance at the Mie frequency lom with a total 
linewidtbJ 5 . 

7t(«) = 7i +7(o)- C 1 ) 

Thus, in addition to the intrinsic linewidth 7;, we have 
to consider a size-dependent contribution which can be 
expressed a a 11 ! 12 ' 25 

** - E« {£) < 2 » 

where e-p = h 2 k F /2m c and vp are the Fermi energy and 
velocity, respectively, and go is a smooth function that 
will be given in (|36|) . 26 The size-dependent linewidth 
7(a) results from the decay of the surface plasmon into 
particle-hole pairs by a Landau damping mechanism, 
which is the dominant decay channel for nanoparticle 
sizes 0.5 nm < a < 5nm considered in this work. For 
larger clusters, the interaction of the surface plasmon 
with the external electromagnetic field becomes the pre- 
ponderant source of damping^ 

In a quantum many-body approach, the surface plas- 
mon appears as a collective excitation of the electron 
system. Discrete-matrix random phase approximation 
(RPA) provides a useful representation since the eigen- 
states of the correlated electron system are expressed 
as superpositions of particle-hole states built from the 
Hartree-Fock ground state. Following similar approaches 
developed for the study of giant resonances in nuclei, 
Yannouleas and Broglia 12 proposed a partition of the 
many-body RPA Hilbert space into a low-energy sector 
(the restricted subspace), containing particle-hole exci- 
tations with low energy, and a high-energy sector (the 
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FIG. 1: Absorption cross section a(ui) in arbitrary units ex- 
tracted from TDLDA calculations for a sodium cluster con- 
taining N = 832 valence electrons. The classical Mie fre- 
quency is um, while com is the frequency of the surface plas- 
mon resonance taking into account the spill-out effect. It has 
been obtained with the Kohn-Sham self-consistent ground- 
state density as explained in Sec. II V CI ^ sp corresponds to 
the position of the maximum of the absorption curve. 



additional subspace). The surface plasmon arises from a 
coherent superposition of a large number of basis states of 
the restricted subspace. Its energy lies in the high-energy 
sector, and therefore the mixture with particle-hole states 
of the additional subspace results in the broadening of the 
collective resonance. 

An alternative to the previous approaches is given 
by numerical calculations using the time-dependent lo- 
cal density approximation (TDLDA) within a jcllium 
modelii2*2i The absorption cross section 

= ^ J! K/N°>! 2 - E f + E o) ( 3 ) 

can be obtained from the dipole matrix elements between 
the ground state |0) and the excited states |/) of the 
electron system with energies Eq and Ef, respectively, 
c is the speed of light. In the absorption cross section, 
the surface plasmon appears as a broad resonance (see 
Fig. ^) centered at a frequency tu sp close to lom and with 
a linewidth approximately described by Q . 

Finally, the separation in center-of-mass and relative 
coordinates for the electron systemSSiSS within a mean- 
field approach allows to describe the surface plasmon as 
the oscillation of a collective coordinate, which is damped 
by the interaction with an environment constituted by a 
large number of electronic degrees of freedom. Within 
such a decomposition, the effects of finite size, a dielectric 
material^ or the finite temperature of the electron gas 
are readily incorporated. Since the excitation by a laser 
field in the optical range only couples to the electronic 
center of mass, this approach is particularly useful. 

The experimentally observed resonance frequencies are 
smaller than the Mie frequency umA^ and this is quali- 



tatively captured by TDLDA calculations since uj sp < ujm 
(see Fig. QJ. This redshift is usually attributed to the 
so-called spill-out effect,^ The origin of this quantum ef- 
fect is a nonzero probability to find electrons outside the 
nanoparticle, which results in a reduction of the effective 
frequency for the center-of-mass coordinate. If a fraction 
N ou t/N of the electrons is outside the geometrical bound- 
aries of the nanoparticle, we can expect that the electron 
density within the nanoparticle is reduced accordingly 
and the frequency of the surface plasmon is given by 



w M = WMyl — (4) 

We will see in Sec.[n]that such an estimation can be for- 
malized considering the form of the confining potential. 
It has been known for a long time^ that u>m > <^sp> and 
therefore the spill-out effect is not sufficient to explain the 
redshift of the surface plasmon frequency as it can be seen 
in Fig. ^ Recently, the coupling to the electronic envi- 
ronment has been invoked as an additional source of fre- 
quency shifti^SiSI In this work we provide an estimation 
of such a contribution (analogous to the Lamb shift of 
atomic physical) and its parametric dependence on the 
particle size and electron temperature (see Sec. 01 . In 
our approach, we always adopt a spherical jellium model 
for the ionic background which is expected to yield reli- 
able results for not too small nanoparticle sizes. Despite 
the fact that we have to make several approximations, 
we find that this additional shift implies a reduction of 
the surface plasmon frequency, in qualitative agreement 
with the TDLDA calculations. However it is known^ that 
the experimental resonance frequency is even lower than 
the TDLDA prediction. This could be due to the ionic 
degrees of freedom. 

In femtosecond time-resolved pump-probe experiments 
on metallic nanoparticles^S the surface plasmon plays 
a key role. The pump field heats the electron system 
via an excitation of the collective mode. The probe then 
tests the absorption spectrum of the hot electron gas in 
the nanoparticle, and therefore it is important to under- 
stand the temperature dependences of the surface plas- 
mon linewidth and of the frequency of the resonance. 

The paper is organized as follows: In Sec.[n] we present 
our model and the separation of the electronic degrees of 
freedom into the center of mass and the relative coor- 
dinates. The center-of-mass coordinate provides a natu- 
ral description of the spatial collective oscillations of the 
electronic cloud around its equilibrium position resulting 
from a laser excitation. Its coupling to the relative coor- 
dinates, described within the mean-field approximation, 
is responsible for the broadening of the surface plasmon 
resonance that we evaluate in Sec. IIIII We determine 
the dependence of the surface plasmon linewidth on the 
size of the nanoparticle and explore its low-temperature 
properties. In Sec. II VI we evaluate the spill-out effect, fo- 
cusing on its dependence on size and temperature of the 
nanoparticle. We show by means of TDLDA calculations 
that the spill-out effect is not sufficient to describe the 
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redshift of the surface plasmon frequency as compared to 
its classical value. In Sec. [V] we propose an estimation 
of the environment-induced redshift of the surface plas- 
mon resonance which adds to the spill-out effect. Both 
effects together, spill-out and frequency shift due to the 
electronic environment, could explain the observed red- 
shift of the resonance frequency. In Sec. IVII we draw 
the consequences of our findings for pump-probe exper- 
iments. The temperature dependences of the linewidth 
and frequency of the surface plasmon resonance peak per- 
mit to qualitatively explain the time dependence of the 
measured optical transmission as a function of the de- 
lay between the pump and the probe laser field. We 
finally conclude and draw the perspectives of this work 
in Sec. EH 



II. COUPLING OF THE SURFACE PLASMON 
TO ELECTRON-HOLE EXCITATIONS 

Treating the ionic background of the nanoparticle as 
a jellium sphere of radius a with sharp boundaries, the 
Hamiltonian for the valence electrons is given by 



N 



H = J2 



i=l 



2m„ 



N 

E 



i 



2 |r, ; 

I J — 1 

(i*7) 



(5) 



where is the position of the i particle and = |r 
The single-particle confining potential 



U(r) 



Ne 2 
2^" 



- 2 - 3a 2 ) Q(a 



9(r - a) (6) 



is harmonic with frequency wm = y/ Ne 2 /m e a 3 inside the 
nanoparticle and Coulomb-like outside. 8(x) denotes the 
Heaviside step function. In principle, the photoabsorp- 
tion cross section (0 can be determined from the knowl- 
edge of the eigenstates of H . However, except for clusters 
containing only few atoms, this procedure is exceedingly 
difficult, and one has to treat this problem using suitable 
approximation schemes. 
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(8) 



is the Hamiltonian for the relative-coordinate system. 

Assuming that the displacement R of the center of 
mass is small compared to the size of the nanoparticle, 
we can expand the last term on the r.h.s. of Q). To 
second order, we obtain 

C/(|r'+R|) -£/(/) ~ R-VC/(r') + i (R • V) 2 U(r'), (9) 

where the derivatives are taken at r = r' (R = 0). Choos- 
ing the oscillation axis of the center of mass in the z- 
direction, R = Ze z , we obtain with © 



R • VU{r') = ZmeUJM 



z'Q{a-r') + — T Q{r' -a) 



and 



(R ■ V) 2 U(r') 



Z 2 Ne 2 



(10) 



(11) 



1 r^, n 1 - 3 cos 2 0' . , ' 
-^&{a - r') H = 9(r' - a) 



Eq. 111)1 represents a linear coupling in Z between the 
center-of-mass and the relative-coordinate system. In 
second order, the first term on the r.h.s. of (|llf> is the 
dominant contribution to the confinement of the center 
of mass. The second term on the r.h.s. of ljll|l is of second 
order in the coupling and therefore is neglected compared 
to the first-order coupling of (|l(Jfl . 

Inserting 10 and (|llfl into Q, we obtain 



H 



lNe< 



N 



2Nm c 



r 2 E ( 



H Icl + H c , (12) 



where 



N 



ff c = ^R-[VC7(K) 



R=0 



(13) 



Separation into collective and relative 
coordinates 



A particularly useful decomposition 28 of the Hamilto- 
nian J5j can be achieved by introducing the coordinate of 
the electronic center of mass R = J2i V i/N and its con- 
jugated momentum P = ^ i p^. The relative coordinates 
are denoted by v\ — — R and = Pi — P/N. Then, 
the Hamiltonian (0 can be written as 



H = 



2Nm 



N 

E 



[£/(|r; + R|)-[/(rO], (7) 



is the coupling between the center-of-mass and the rela- 
tive coordinates to first order in the displacement R of 
the center of mass according to (1101) . The remaining sum 
over i in (|12|) yields the number of electrons inside the 
nanoparticle, i.e., N — N ont where JV ou t is the number 
of spill-out electrons, and finally we rewrite the Hamilto- 
nian as 

H = H cm + H le} + H c . (14) 
The Hamiltonian of the center-of-mass system is 

p2 j 



2Nm c 



(15) 
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It is the Coulomb tail of the single-particle confinement 
© which yields the frequency wm given in Q instead of 
lum for the effective harmonic trap that is experienced by 
the center-of-mass system. 

Eq.[n]with ©, {E3), and {El recovers up to the sec- 
ond order in R the decomposition derived in Ref. 1281 
However, in contrast to that work, we do not need to ap- 
peal to an effective potential for the center-of-mass sys- 
tem. 

The structure of (|I4|I is typical for quantum dissipative 
systems s24 The system under study (H cm ) is coupled via 
H c to an environment or "heat bath" H IC \, resulting in 
dissipation and decoherence of the collective excitation. 
In our case the environment is peculiar in the sense that 
it is not external to the nanoparticle, but it represents a 
finite number of degrees of freedom of the gas of conduc- 
tion electrons. If the single-particle confining potential U 
of © were harmonic for all r, Kohn's theorem^ would 
imply that the center of mass and the relative coordinates 
are decoupled, i.e., H c = 0. Thus, for a harmonic poten- 
tial U the surface plasmon has an infinite lifetime. The 
Coulomb part of U in 10 leads to the coupling of the cen- 
ter of mass and the relative coordinates, and translates 
into the decay of the surface plasmon. Furthermore, it 
reduces the frequency of the center-of-mass system from 
wm according to Q. 



B. Mean- field approximation 



Introducing the usual annihilation operator 
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FIG. 2: LDA self-consistent potential V as a function of 
the radial coordinate r for a sodium nanoparticle containing 
N = 1760 valence electrons. The radius a is indicated by the 
vertical dotted line. The Fermi level is marked by the dashed 
line. 



r — a. However, the details of the self-consistent poten- 
tial close to the surface may be crucial for some proper- 
ties, as we show in Sec. IIV CI 

Inserting l|l(J|l into the coupling Hamiltonian H c i|lo|l . 
and expressing the .^-coordinate in terms of the creation 
and annihilation operators of l|16l) , one obtains in second 
quantization 



a/3 



(19) 



Nm c ujM 



2h 



y/2Nm c hi~L!M 



P? 



(16) 



where 



where Pz is the momentum conjugated to Z, and the 
corresponding creation operator w, H cm reads 



H cm = huj M tfb 



(17) 



Since the electron-electron interaction appears in the 
Hamiltonian JHJl, it is useful to describe H TC i within a 
mean-field approximation. One can write 



(18) 



where e a are the eigenenergies in the effective mean-field 
potential V and cj, (c a ) are the creation (annihilation) 
operators associated with the corresponding one-body 
eigenstates \a). The mean- field potential V can be de- 
termined with the help of Kohn-Sham LDA numerical 
calculations. In Fig. [21 we show it as a function of r/ao 
where do ~ 0.53 A is the Bohr radius, for a nanoparticle 
containing N = 1760 valence electrons. We see that it 
is relatively flat at the interior of the nanoparticle, and 
presents a steep increase at the boundary. The potential 
jump is often approximated by a true discontinuity at 



d a R = (a | 



z@(a — r) H r 9(r - a) 



\/3) (20) 



is the matrix element between two eigenstates of the un- 
perturbed mean- field problem. In (|19fl . we have defined 
the constant A = y l /hm c ujy l /2N and neglected the spill- 
out for the calculation of the coupling when we expressed 
Z in terms of b and w . The sums over a and (3 in (|19fl are 
restricted to the additional (high energy) RPA subspace 
mentioned in the introduction. If initially the center of 
mass is in its first excited state, the coupling H c allows 
for the decay of the collective excitation into particle-hole 
pairs in the electronic environment -ffrci, the so-called 
Landau dampingii&SiSi 



III. TEMPERATURE DEPENDENCE OF THE 
SURFACE PLASMON LINEWIDTH 

In this section, we address a semiclassical calculation 
of the surface plasmon linewidth and extend our zero- 
temperature results of Ref. to the case of finite tem- 
peratures. We are interested in nanoparticles with a large 
number of confined electrons, typically of the order of 10 2 
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to 10 4 . Thus, the Fermi energy is much larger than the 
mean one-body level spacing and the semiclassical ap- 
proximation can be applied. 



A. Fermi's golden rule 

When the displacement of the center of mass is much 
smaller than the size of the nanoparticlc, it is possible 
to linearize the coupling Hamiltonian as it was done in 
Sec.m The weak coupling regime allows to treat H c as a 
perturbation to the uncoupled Hamiltonian H cm + H re \. 
Assuming that initially the center of mass is in its first 
excited state |l cm ), i-e., a surface plasmon is excited, two 
processes limit the surface plasmon lifetime: the decay 
into the ground state with creation of a particle-hole pair 
and the excitation to the second state of the center of 
mass accompanied by the annihilation of a particle-hole 
pair. Obviously, this last process is only possible at fi- 
nite temperatures. Then, Fermi's golden rule yields the 
linewidth 

7= 1T E ^.K^cm^rell-ffcllcm^rel)! 2 
/,ci,F,d 

X 5(M> M - £ Fcm + £/ rol ~ SFrel) ( 21 ) 

for the collective state |l cm )- ln the golden rule, \F cm ) 
and l-Frei) are the final states of center-of-mass and rela- 
tive coordinates with energy e p cin and e p Tei , respectively. 
The probability of finding the initial state |J re i) occupied 
is given in the grand-canonical ensemble by the matrix 
element 



(/ rcl |e-^--^)|/ rcl ) 



(22) 



of the equilibrium density matrix at the inverse temper- 
ature (3 = l/k^T, fee being the Boltzmann constant, 
/i is the chemical potential of the electrons in the self- 
consistent field V and S the grand-canonical partition 
function. Introducing the expression l|19|) of the coupling 
Hamiltonian, one finds that 



7 = S(o)m) + 2E(-cD M ), 



(23) 



where the first term is the rate associated with the spon- 
taneous decay of a plasmon, while the second one is the 
rate for a plasmon excitation by the thermal environ- 
ment. The relative factor of 2 arises from the dipolar 
transition between the first and second center-of-mass 
excited states, (2 cm |6 1 '|l cm ) = 

In i|23|) . we have introduced the function 



plasmon induced by the electronic environment will be 
determined. In the above expression, tOap = (e a — £p)/h 
is the difference of the eigenenergies in the self-consistent 
field V and 



i 



(25) 



is the Fermi-Dirac distribution. It is understood that 
in (|24fl . \uj a p | > uj c where u c is some cutoff separating 
the restricted subspace that builds the coherent super- 
position of the surface plasmon excitation from the ad- 
ditional subspace at high energies. The expression Ij24(l 
implies the detailed-balance relation 



E(-w) = e-^EH 
which allows to write (|23|) as 



7 



= E(w M ) (1 + 2c 



(26) 



(27) 



This expression shows that an excitation of the surface 
plasmon to a higher level is suppressed at low tempera- 
tures. 



B. Semiclassical low-temperature expansion 

We now use semiclassical techniques to determine 
the low-temperature behavior of E(w). In view of the 
detailed-balance relation Ij26(l we can restrict ourselves 
to positive frequencies lu. In a first step, we need 
to calculate the matrix elements d a p defined in l|2(J|l . 
The spherical symmetry of the problem allows to write 
d a p = A] 7 ^ rn ' 3 TZ(e a ,ep), where an expression for the an- 
gular part •A™"!™' 3 is given in Ref. |2^ Dipole selection 
rules imply l a = lp ± 1 and m a = mp, where I and m 
are the angular momentum quantum numbers. The self- 
consistent potential is usually approximated by a step- 
like function, V(r) ~ Vo0(r — a), where Vb = £f + W 
with W the work function of the metal. The accuracy of 
such an approximation can be estimated by comparison 
with LDA numerical calculations (see Fig.|2and Ref.|29J). 
In the case of strong electronic confinement Vq 3> £p, the 
radial part can be approximated byi£ 



TZ(e a ,ep) 



2h 2 



m c a ( Ea - erf 



(28) 



The condition of strong confinement implicit in l)28|) as- 
sumes that the spill-out effect is negligible. Replacing in 
(|24|l the sums over a and [3 by integrals and introduc- 
ing the angular-momentum restricted density of states 
(DOS) ei (e)^&.one gets 



47T f°° 



f(e)]f{e-Hu) 



(29) 



a/3 

(24) 

which will be helpful for the evaluation of 7 and will 
also be useful in Sec. Ivl when the redshift of the surface 
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Qi(e)gi>(e-tiw) AJfiT R(e,e - fiu) 



6 



where a factor of 2 accounts for the spin degeneracy. 

We now appeal to the semiclassical approximation for 
the two DOS appearing in (|29(l using the Gutzwiller trace 
formula^ for the effective radial motioniSSi The DOS is 
decomposed into a smooth and an oscillating part. With 
the smooth part 



yf2m c a?e/h 2 -{l + l/2f 
2tt£ 



(30) 



of the DOS, and after performing the summation over 
the angular momentum quantum numbers, (|29[l reads 



E(w) 



\ 3 

3«p / \ £f 



8a \ £f 
Here, we have introduced 



^M = /JgH-/( e )]/( e -^(^ 
where 

H(x) = (2x - l)y/x{x- 1) - In (y/x + y/x~T) 



(31) 



(32) 



(33) 



is an increasing function. The dependence on the chemi- 
cal potential /i in i|32|) is via the Fermi functions appear- 
ing in this expression. We calculate the function T in the 
appendix and this leads for temperatures much smaller 
than the Fermi temperature Tp to 



£f J_ 



with 



.'/ ( r, — ) = //,,(•'•) -f fi- 2 U) 



T 




= ^{y/x(x + i)(4x(x + l) + 3) 

- 3(2x + 1) In (y/x~ + Vx + 1) 

+ 3) 



(34) 



(35) 



(36) 



yj x{x — 1) (Ax{x — 
- Z{2x - 1) In (y/x + y/x^T) 0(x-l)} 

is a monotonically increasing function with go(0) = and 
lim^oo g (x) = 1, while 

7r 2 r 

ff2 (a:) = —{y/x(x + l)(6x - 1) + In (y/x + Vx + l) 

- \/xJx - l)(6x + 1) + In (y/x + y/x-l) 9(x-l)| 

(37) 

with (72(0) = and lima, _»oo g^ix) — ir 2 /6. For x near 1, 
a nonanalytical T 5 / 2 -correction coming from i|A5|l has to 
be added to 



Eq.E3then yields the size- and temperature-dependent 
surface plasmon linewidth 



7 



3^F 

la"' 



e F T 
hujyi ' Tf 



(38) 



Note that the exponential factor appearing in (|27|l is ir- 
relevant in our low-temperature expansion. Furthermore, 
we have replaced (Dm by lum for the calculation of 7. As 
we will see in Sec. lIVI the spill-out correction to the Mie 
frequency scales as 1/a. Furthermore, (|38|l shows that 
7 scales also as 1/a. Thus, incorporating the spill-out 
correction in the result of l|38|l would yield higher order 
terms in 1/a, inconsistently with our semiclassical expan- 
sion that we have restricted to the leading order. 

At T — 0, we recover with l)38|) the well-known 1/a 
size dependence of the surface plasmon linewidth. First 
found by Kawabata and Kubo^ this size dependence 
is due to the confinement of the single-particle states in 
the nanoparticle. We also recover the zero-temperature 
frequency dependence found in Refs. 11 and 12; As a 
function of ep, the Landau damping linewidth increases 
linearly for £p <C hu>M and as ^/Fp for £p ^> Hlum- The 
increase of the linewidth is a consequence of the fact 
that with increasing Fermi energy the number of rele- 
vant particle-hole excitations rises. 

Since the function g2 is positive, finite temperatures 
lead to a broadening of the surface plasmon resonance 
which to leading order is quadratic. As for T = 0, the 
linewidth decreases with increasing size of the nanopar- 
ticle like 1/a. Based on classical considerations, this has 
been proposed in Ref. Il2l where the authors argued that 
7 ~ v/a with the average speed 



3«F 



1 



2 



(39) 



of electrons at the temperature T. The result of Ref. [12J 
is only relevant in the high energy limit ep S> S^m where 
it agrees with our general result l|38|). 

In Fig. we show the linewidth (|38|l as a func- 
tion of the temperature, scaled by the zero-temperature 
linewidth, for different values of the ratio £f/Swm (solid 
lines). In order to confirm the validity of our low tem- 
perature expansion, we compare it to the result of a nu- 
merical integration (dashed lines) of (|27|) with Eqs. 1311 - 
13*31 The agreement of the expansion (|38|l with this direct 
integration is excellent for low temperatures. For fixed 
T/Tp, the deviation increases with £ F /Swm since then, 
the low-temperature condition Hlum S> k^T is less and 
less fulfilled. 

There are a number of experiments showing a 
broadening of the surface plasmon linewidth with the 
temperatureiSSiiil Those experimental results on not 
too small clusters indicate the presence of small correc- 
tions to the width of the collective surface plasmon exci- 
tation due to finite temperatures, in agreement with our 
result (|38|l . In Ref. 38, absorption measurements on gold 
nanoparticles with a diameter ranging from 9 to 25 nm in 
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FIG. 3: Surface plasmon linewidth as a function of the tem- 
perature, scaled by the linewidth at T = 0, for different ra- 
tios £F/ftajM- Solid lines: low-temperature expansion 1381 . 
For £f = huoyi, we have accounted for the additional y 5 / 2 _ 
correction to 7, see l)A5)l . Dashed lines: numerical integration 
of 12711 with Eqs. 1311)331 The arrow indicates the direction of 
increasing £f/^m. 



aqueous solution have shown only a weak temperature ef- 
fect on the surface plasmon linewidth. In Ref. 39, a weak 
broadening of the plasmon resonance in silver and gold 
nanoparticles of sizes a = 4.25 to 10 nm is reported, ac- 
companied by a small redshift of the peak position as the 
temperature increases. In Ref. lifj the temperature de- 
pendence of small silver clusters (radius of 1.6 to 10.5 nm) 
embedded in a glass matrix has been investigated and a 
rather small broadening of the plasmon line has been re- 
ported as the temperature increases from 1.5 to 300 K. 

There are also experiments on very small clusters^ 
showing a strong temperature effect on the surface plas- 
mon linewidth. It has been observed in Ref. |6J that the 
plasmon width of mercury clusters increases dramati- 
cally with temperature. A systematic study of the tem- 
perature dependence of the linewidth in small charged 
sodium nanoparticles with N = 8, 20, and 40 valence 
electrons has been carried out in Ref. 0. As the tempera- 
ture of the cluster is increasing, the authors found a pro- 
nounced broadening of the resonance which goes typically 
as Vt. This is in apparent contradiction with our result 
(J3SJ. However, for very small particle sizes, an additional 
broadening mechanism becomes important, namely the 
coupling of the surface plasmon to quadrupole surface 
thermal fluctuations 

In addition to the smooth part 1)30)1 of the semiclassical 
DOS, there exists a term which oscillates as a function of 
the energy. The convolution of the two oscillating parts of 
the DOS in 1)29)1 yields an additional contribution to the 
plasmon linewidth, which oscillates as a function of the 
size of the nanoparticle. This contribution is significant 
only for very small sizes. Its existence has been con- 
firmed by TDLDA calculations. 29,36 Size-dependent os- 
cillations are also observed in experiments on very small 
nanoparticlesii 3 . The oscillations of the linewidth can 



be expected to be smoothed out with increasing tem- 
perature because of thermal broadening suppressing the 
particle and hole oscillations of the DOS. 



IV. SPILL-OUT-INDUCED FREQUENCY SHIFT 

We now turn to the evaluation of the redshift of the 
surface plasmon frequency with respect to the classical 
Mie value wm- In this section, we calculate one of the 
contributions affecting the resonance frequency, namely 
the spill-out effect, while in Sec. we will calculate the 
additional frequency shift induced by the electronic en- 
vironment. 



A. Mean-field approximation 

Obtaining the many-body wave function and extract- 
ing the number of spill-out electrons is very difficult al- 
ready for very small clusters, and practically impossible 
for larger particles. Therefore we use the mean-field ap- 
proximation and treat the electronic degrees of freedom 
in the mean-field one-particle potential V shown in Fig. [21 
Within this approximation, the number of electrons out- 
side of the nanoparticle is given by 



AUt = 



2/ d E y>(e)/( £ ) / d 3 r|^ ;m (r)| 

JO lm J(r>a) 



(40) 

where Qi{e) is the DOS restricted to a fixed angular mo- 
mentum I from ()30)l . The factor of 2 accounts for the 
spin degeneracy. Because of the spherical symmetry of 
the problem, the one-particle wave function 



(41) 



separates into a radial and an angular part given by the 
spherical harmonics Yj m {9, ip). The radial wave functions 
u £ i(r) satisfy the reduced Schrodinger equation 



h 2 d 2 h 2 m 

+ 



1 



2m„ dr 2 



2m r r 2 



V(r) 



u £ i{r) = eu e i{r) (42) 



with the conditions u e ;(0) = and linv^oo [u e i{r)/r] — 
0. Eq. 021 yields the single-particle eigenenergies e in the 
mean-field potential V, that we approximate by the step- 
like potential V(r) = Vo&{r — a). 

The Fermi function in 1)40) suppresses contributions to 
the energy integral from values higher than £f plus a few 
k-&T. For low temperatures k-gT <C W, the states in the 
continuum do not contribute to N ou t and we restrict our 
evaluation to the bound states with e < Vb- 44 Defining 
k 



y/2m c s/ h and x = \/2m e (Vb 
regular solutions of l|4*2")l 



s)/H, we find for the 



u £ i{r) 



r j A k iVkJ l+ i(kr), r ^ a, 
Bkiy/xK l+ i{xr), r>a, 



(43) 
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where J v (z) are Bessel functions of the first kind and 
K v (z) are modified Bessel functions. The normalization 
constants Am and Bm are given by 



with 



Cm 



A m = 
Bm = 

^+i( fcQ ) 
K l+ i{xa) 



kaCh 



Ji+d ka ) 



XaC k i K l+ i(xa)' 



(44a) 
(44b) 



^_i(xa)^+|(xa) 
- Ji_x{ka)J l+ 3(ka). 



(45) 



We therefore obtain for the integral in 140f) 



(r>o) 



d 3 r|^ /m (r)|' 
Jf +X (ka) 



a 



M 



- 1 



(46) 



The summation of this expression over all one-particle 
states required to obtain N out according to H4U|) cannot 
be done exactly. We therefore use a semiclassical ap- 
proximation which provides additional physical insight 
into the spill-out effect. 



B. Semiclassical low-temperature expansion of the 
number of spill-out electrons 

The integral of the electronic density (|46|) increases 
with the energy. Combined with the increasing DOS in 
(|40|l and the Fermi function providing an energy cutoff, 
this allows to conclude that the spill-out is dominated by 
the energies near the Fermi level. In addition, in most 
of the metallic nanoparticles, we have £p ~ W ^> A, 
where A is the mean single-particle level spacing. The 
semiclassical limit, in which ka and x a must be much 
larger than one, then applies. In this limit, we obtain for 
the integral l|46|) of the density outside of the nanoparticle 



d 3 r \lpelm( r ) 



(r>a) 



XaV 



(47) 



to first order in 1/ka and l/x a - Since this result does not 
depend on the angular momentum quantum numbers I 
and m, the total DOS V, pi(e) (whose explicit expres- 
sion can be found in Ref. |29|) is sufficient to determine 
the number of spill-out electrons. Inserting (|47(l into 1(40(1 
then yields 



Am e a 
3^¥Vq 



v 



def(s) 



r 3/2 



(48) 



Applying the low-temperature Sommerfeld expansion, 45 
we obtain 



with 

where 

Co (a:) 
and 



_ (fcpfl) fep T_ 

67T \V T F 



l r 



-y/x(l-x)(2x + 3)+3 



arcsmvi 



2 / \ 3/2 

7T / X 



3 \l-x 



(2-1). 



(49) 



(50) 



(51) 



(52) 



Note that the upper bound of the integral over the en- 
ergy in 1(4*51) has been replaced by Vq since we neglect 
exponentially suppressed contributions from higher ener- 
gies. Therefore our Sommerfeld expansion is reliable for 
temperatures T/Tp < 1 — Sf/Vq. 

At zero temperature the number of spill-out electrons 
increases smoothly with increasing Fermi energy, reach- 
ing its maximal value (fcpa) 2 /4 at Ep =Vq. According to 
(|49|l and H5U|) , the number of spill-out electrons increases 
with temperature. This is expected since the evanescent 
part of the wave function increases with the energy of the 
occupied states. 

Scaling the result 149(1 with the total number of elec- 
trons in the nanoparticle, N — 4(fcpct) 3 /97r, we obtain 



N 



8fc F a \VqT ¥ ) 



(53) 



This relative spill-out scales to first order as 1/a, and 
is therefore negligible for large particles. As the work 
function W depends on the size of the nanoparticle like 
W = + a/ar^ the number of spill-out electrons i(53() 
acquires, via the dependence on Vo, higher-order correc- 
tions to the 1/a scaling. We neglect these terms because 
they are of the same order as terms neglected in our semi- 
classical expansion. Therefore we approximate W by its 
bulk value Woo- The dependence of (|53|l on a can be 
interpreted by observing that the spill-out is a surface 
effect so that 7V ou t increases only with a 2 . Inserting llo.'il) 
into (0J, one can calculate the spill-out-induced redshift 
of the surface plasmon resonance. This redshift increases 
for decreasing sizes and for increasing temperatures, in 
qualitative agreement with experimentsi^2ii£ 
One can define a spill-out length as the depth 



, 1 ^Vout 

ls= 3— a 



(54) 



of the spill-out layer. Inserting 1(53(1 into 1(54(1 yields the 
size-independent spill-out length 



(55) 
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FIG. 4: k F l s from 1551 1 as a function of ef/Vo at zero tem- 
perature (black line) and for finite temperatures (solid gray 
lines). The dashed lines result from the numerical integration 
of 14811 . In the figure, T/T F increases from bottom to top. 



which is represented in Fig. 0] as a function of the ratio 
£f/Vo for different temperatures. The result is compared 
with the results of a numerical integration of (|48|l (dashed 
lines). This confirms our expectations about the valid- 
ity of our result, namely for low temperatures and for a 
reasonable ratio £f/Vq. 



Number of spill-out electrons: 
vs. LDA 



Semiclassics 



In this section, we compare our semiclassical evalua- 
tion of the spill-out effect at zero temperature with LDA 
calculations on spherical sodium nanoparticles. One pos- 
sible way to estimate -/V out is to use the LDA, which al- 
lows to compute the spherically symmetric self-consistent 
electronic ground-state density for nanoparticles with 
closed electronic shells. Integrating the density outside 
the nanoparticle then yields an approximation to A out , 
and thus to l s according to l|54|). An estimation of the 
spill-out length from (|55|l by means of our semiclassical 
theory at zero temperature gives l s ~ 0.2 an for sodium 
nanoparticles, while Madjet and collaborators obtained 
on the basis of Kohn-Sham and Hartree-Fock calculations 
l s around 0.55 an for clusters of size N = 8-196^ With 
our LDA calculations, we obtain l s of the order of 0.45 an 
for all closed-shell sizes between N — 8 and N = 1760 
(see the squares in Fig. 0. 

The fact that the semiclassical spill-out length is sig- 
nificantly smaller than that of LDA is a consequence of 
our assumption of a step-like potential for V. Indeed, the 
LDA self-consistent potential V shown in Fig. [5] deviates 
from the form V(r) = Vo®(r — a) that we have used. As 
one can see in Fig. [SJ the Fermi level does not coincide 
with V(a). Defining the effective radius of the nanopar- 
ticle for the spill-out effect by means of V(a c ff) = £p, 
it seems appropriate to approximate the self-consistent 
potential by V(r) = Vo0(r — a e s). An estimation from 
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FIG. 5: Spill-out length l E at T = in units of the Bohr 
radius ao as a function of the total number of electrons N, 
after 15511 (solid line) for sodium nanoparticles. The dashed 
line is twice the result of 155H . and is obtained considering 
the effective radius a e a for the approximated self-consistent 
potential, namely V(r) = Vo@(r — a c g). The squares result 
from LDA calculations. 



our LDA calculations gives a G ff ~ a + l s for all sizes be- 
tween N = 8 and N = 1760. Using this effective radius 
in H49fl does not change our results for AT out and l s since 
a^> l s . However, the spill-out length is defined from the 
geometrical radius a of the ionic jellium background of 
the nanoparticle. Since a c g — a + l s , it actually yields 
an effective spill-out length ^ ff « 2l s and approximately 
doubles our result for -/V out . The improved T — result 
for the spill-out length is represented by the dashed line 
in Fig. |3J It yields good agreement with the spill-out 
length as deduced from the LDA calculations (squares in 
Fig.EJ. Note that there is no need to consider an effective 
radius of the nanoparticle in the calculation of the reso- 
nance width presented in Sec. IIIII Indeed, the linewidth 
7 scales as 1/a. Replacing a by an effective radius for the 
linewidth 7 would thus only lead to higher-order correc- 
tions. 



D. Redshift of the surface plasmon resonance 

We now examine the redshift of the surface plasmon 
resonance by means of TDLDA calculations. In Fig. [7| 
we show the frequencies deduced from the LDA number 
of spill-out electrons for various closed-shell nanoparti- 
cle sizes between N = 8 and N — 1760, where N out 
is incorporated according to (squares). The dashed 
line is our semiclassical result from (I53|) with J3J, where 
we have taken into account the effective radius for the 
self-consistent potential, as discussed in the preceding 
section. We see that our analytical expression for the 
spill-out is in a good agreement with the LDA calcula- 
tions, and that the redshift is increasing with decreasing 
size as predicted by 

Alternatively the resonance frequency can be extracted 
directly from the absorption cross section cr(uj) defined in 
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FIG. 6: Detail of the LDA self-consistent potential V close 
to the boundary of the cluster as a function of the radial 
coordinate r for a sodium nanoparticle containing N = 1760 
valence electrons. The radius a is indicated by the vertical 
dotted line, and the effective radius defined by V{a e s — a + 
l s ) = £f is indicated by the dashed line. The Fermi level 
corresponds to the horizontal dashed line. 




FIG. 7: Surface plasmon frequency as a function of the 
number N of valence electrons in the nanoparticle. Squares: 
iIim/um, frequencies deduced from the LDA number of spill- 
out electrons according to Q. Dashed line: semiclassical eval- 
uation of 153t with Dots: ui sp /lum, frequencies obtained 
by fitting the TDLDA absorption curves with a Mie-type cross 
section?" Dotted line: linear fit to the dots. Solid line: sum of 
the spill-out effect (dashed line) and the environment-induced 
redshift of the resonance (16411 . All results shown here are at 
zero temperature for the case of sodium. The points shown 
represent all closed shell sizes between N = 8 and N = 1760. 



© and calculated from the TDLDA response function 
(see Fig.nj. An upper bound for the resonance energy is 
given by the root mean squar o 2 ! 47 



r— J dujuj 2 a{lS) 

v ( w ) = \ -ran — 
V Jo duj 



(56) 



plasmon frequency from the Mie value. In the spherical 
jcllium model, the frequency deduced from l|5()[) coincides 
with CuMr^ the Mie frequency redshifted by the spill-out 
effect (0J. We have checked with our TDLDA calcula- 
tions that the frequencies deduced from N ont , i.e., from 
the electronic ground-state self-consistent density corre- 
spond to the frequencies obtained from l|56|) . up to the 
numerical error. At zero temperature, the TDLDA thus 
fulfills the sum rule^ -J (to 2 ) = Com- 

The dots in Fig. \7\ represent the maxima tu sp of the 
absorption curves obtained by fitting the TDLDA results 
with a Mie-type cross sectioni It is clear from Fig. [7| 
that the frequency obtained by J3J from the number of 
spill-out electrons (squares) overestimates the resonance 
frequency u> sp (see also Fig. ^| It has been noticed in 
Ref. that the spill-out effect is not sufficient to describe 
the experimentally observed resonance frequency. 

In the following section, we propose to interpret the 
discrepancy between w sp given by the TDLDA and Com 
deduced from the spill-out by means of the coupling of 
the surface plasmon mode to electron-hole excitations. 
This coupling results in a shift of the surface plasmon 
frequency which adds to the effect of the spill-out. 



ENVIRONMENT-INDUCED FREQUENCY 
SHIFT 



In the absence of the coupling H c , the energy of the 
n th eigenstate of the center-of-mass system is given by the 
eigenenergies £n^ — uHCom of H cm given in ljT7J) . How- 
ever, the coupling Hamiltonian (|19f) perturbs the eigen- 
states of H cul . The leading contribution to the resulting 
shift of the eigenenergies £n is determined using per- 
turbation theory in H c . While the first order does not 
contribute due to the selection rules contained in the cou- 
pling (|19fl , there are four second-order processes involving 
virtual particle- hole pairs. In addition to the two reso- 
nant processes mentioned in Sec. IIII Al i.e., decay into 
the ground state with creation of a particle-hole pair and 
excitation to a higher state accompanied by the annihi- 
lation of a particle-hole pair, there exist also two antires- 
onant processes. A plasmon can be excited to a higher 
collective state by creating a particle-hole pair, or a plas- 
mon can decay by destroying a particle-hole pair. Taking 
into account all four processes, we obtain the resonance 
energy to second order in the coupling 



hid* 



(2) 



(2) 



h {Cj m - 5) 



(57) 



with the frequency shift 

K^cm; -^rcl \H C | lcim Zrel) 



providing a lower bound for the redshift of the surface 



I \F cm , F Te \ | H c 1 cm j -fid , 
£ F cm + £F rt ., - 



(58) 
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where we use the same notations as in Sec. IIII Al Writ- 
ing explicitly the coupling Hamiltonian H c (|19|) in the 
preceding expression, we obtain 



S = /fas) IM»,j| 9 ■ 

a/3 



fi 2 ' 



(59) 



J a/3 



M 



Here, P denotes the Cauchy principal value. The reso- 
nance frequency w sp l|57|l thus contains a correction with 
respect to the value induced by the spill-out effect. As 
we will show, the shift 6 is positive, and thus redshifts 
the plasmon resonance from (Dm- 

According to (|24[) and l|59|l . the energy shift 8 is re- 
lated to the function £(u>) through the Kramers-Kronig 
relation 



1 



+°° , wE(w) 

duj—z 7T 



(60) 



In H6U|) . the frequency cDm appearing in the denomina- 
tor can be replaced by u>m- Indeed, the function £(w) 
of is proportional to l/k-pa. Thus, taking into ac- 
count the spill-out effect in the evaluation of 8 would 
yield higher order terms in powers of 1/fcpa, that we ne- 
glect in the semiclassical limit. Furthermore, we have to 
restrict the integral over the frequency w by introduc- 
ing the cutoff uj c discussed in Sec. IIII Al It arises from 
the fact that the particle-hole pairs that participate to 
5 in 159|) belong to the high-energy sector of the RPA 
Hilbert space, while the surface plasmon excitation is the 
superposition of particle- hole pairs of the restricted low- 
energy subspace. The TDLDA absorption cross section 
shows a large excitation peak at the frequency uj sp which 
supports almost all of the dipole strength. This peak is 
surrounded by particle-hole excitations that become no- 
ticeable for frequencies larger than ~ com — "HI, where rj 
is a constant of the order of unity. Thus for the purpose 
of calculating the integral l|6UI) we can take the cutoff at 
w m — Vl ancl approximate 



(61) 



For frequencies to larger than wm — ?77, the function 
g appearing in l|34|) can be replaced by its asymptotic 
expansion for ep < fiw, 



.9 
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Tf~ 
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2tt2 
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T 



(62) 



Inserting this expression into (|6ip. and performing the 
remaining integral over u> in the semiclassical limit kpa^S> 
1, we arrive at 
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(63) 



We remark that the dependence of this result on the cut- 
off is only logarithmic. 

Inserting our expression H38[l for the linewidth 7, we 
obtain to second order in T/Tp 



where 



3i>f . / £f T 
4a \ fkJM ' Tf 



with 



4v^ 



15tt 



In 



/ 8fcpa 
V 3 Wo (x) 
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h{x) = 



■kJx 



In 
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(64) 



(65) 



(66) 



12 g 2 (x ) 
5?r 2 g (x 



(67) 



The functions go and gi are defined in l|36l) and l|37f) . re- 
spectively. While the linewidth (|38() goes as 1/a, the fre- 
quency shift scales as 1 /a up to a logarithmic factor. This 
redshift increases with temperature and adds to the red- 
shift arising from the spill-out effect discussed in Sec. II VI 

The importance of the shift S can be seen in Fig. 
There, the position of the surface plasmon resonance 
peak for sodium clusters calculated from TDLDA (dots) 
is in qualitative agreement with our semiclassical result 
(solid line) taking into account the shift 8 and the spill- 
out (see Eqs.HandEHlin Sec.|lV|). 50 In Fig.0 we have 
used 77 = 1/2 for the frequency cutoff in l|61|) . Our ap- 
proximate expression for 5 does not allow us to obtain a 
quantitative agreement with the position lo sp of the reso- 
nance frequency shown by the dots. This is not surprising 
considering the approximations needed in order to derive 
an analytical result. First, the expression for 5 is based 
on our result (|34(l for E(w) which was derived under the 
assumption of perfectly confined electronic states. Thus, 
the delocalized self-consistent single-particle states have 
not been accurately treated in the calculation of i5. Sec- 
ond, the cutoff introduced above is a rough estimate of 
the energy beyond which particle-hole excitations couple 
to the surface plasmon. Despite those approximations, 
our estimate implies an increase of the redshift beyond 
that caused by the spill-out (dashed line in Fig. [7J . Com- 
paring the two effects leading to a redshift of the surface 
plasmon frequency, we find that they have the same size 
and temperature dependence, and are of the same or- 
der of magnitude. Therefore, one has to take into ac- 
count both contributions in quantitative descriptions of 
the surface plasmon frequency. 

For zero temperature, the shift 8 has also been consid- 
ered in Refs. I281 andl30l The authors of Ref . l28l have used 
the separation of the collective center-of-mass motion 
from the relative coordinates. Their coupling between 
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the two subsystems which is only nonvanishing outside 
the nanoparticle leads to a shift that they have numeri- 
cally evaluated by means of the RPA plus exc hang e in the 
case of small charged sodium clusters. In Ref.[3(j, the au- 
thors assumed a certain expression for the coupling, and 
a variational RPA calculation was used to obtain an an- 
alytical expression of the environment-induced redshift. 
In contrast to our findings, Refs. [28j and [30J obtained a 
shift 6 proportional to the number of spill-out electrons. 

For very small clusters (N between 8 and 92), a non- 
monotonic behavior of the resonance frequency as a func- 
tion of the size of the nanoparticle has been observed 
experimentally^ This behavior is in qualitative agree- 
ment with our numerical calculations (see dots in Fig. 0) 
and can be understood in the following way. We have 
shown in Refs. [29 and 3^ that the linewidth of the surface 
plasmon resonance presents oscillations as a function of 
the size of the nanoparticle. Furthermore, we have shown 
in this section that the linewidth and the environment- 
induced shift are related through the Kramers-Kronig 
transform ((60(1 . Thus, the shift 5 should also present os- 
cillations as a function of the size of the cluster. This is in 
contrast to the spill-out effect discussed in Sec. llVl where 
the oscillating character of the DOS leads to a vanishing 
contribution, as confirmed by the LDA calculations (see 
squares in Fig.|7|). Those significantly different behaviors 
could permit to distinguish between the two mechanisms 
contributing to the redshift of the surface plasmon fre- 
quency with respect to the classical Mie value. 

Using temperature-dependent TDLDA calculations, 
Hervieux and Bigot^i have recently found a nonmono- 
tonic behavior of the energy shift of the surface plasmon 
frequency as a function of the temperature for a given 
nanoparticle size. They observed a redshift of the surface 
plasmon resonance up to a certain critical temperature 
(e.g., 1000 K and 2500 K for Nai 38 and Na^ 39 , respec- 
tively), followed by a blueshift of the resonance at higher 
temperatures. This behavior is not present in our the- 
ory. The authors of Ref. attribute the nonmonotonic 
temperature dependence to the coupling of the surface 
plasmon to bulklike extended states in the continuum, 
which causes a blueshift, as one can expect for a bulk 
metal. However, we have restricted ourselves to low tem- 
peratures compared to the work function of the nanopar- 
ticle, where we can neglect those extended states in the 
evaluation of the spill-out. The critical temperature of a 
thousand degrees is much smaller than the Fermi temper- 
ature for metals, such that our treatment should remain 
a good approximation. 



VI. TIME EVOLUTION OF THE OPTICAL 
TRANSMISSION IN A PUMP-PROBE 
CONFIGURATION 

We now examine the experimental consequences of 
the temperature dependence of the Landau damping 
linewidth 1(38(1 and of the energy shifts induced by the 



spill-out fEas. HI and 153(1 and by the electronic environ- 
ment 

We focus on the absorption cross section defined in ||SJ . 
Assuming it to be of the Breit-Wigner form, we have 

T) = s(a) 1^11 (68 ) 

>-c sp (T)] 2 + [ 7 (T)/2] 2 ' 

where w sp (T) is the temperature-dependent resonance 
frequency of the surface plasmon excitation given in 157( 1 . 
s(a) is a size-dependent normalization prefactor. For 
the spill-out, we consider the effective radius discussed 
in Sec. IIV(J1 which doubles the results of JSJ). 

Unfortunately, to the best of our knowledge, system- 
atic experimental investigations of the shape of the ab- 
sorption cross section as a function of temperature are 
not available. However, an indirect approach is offered 
by pump-probe experiments, 32 -S where the nanoparti- 
cles are excited by an intense laser pulse. After a given 
time delay, the system is probed with a weak laser field, 
measuring the transmission T of the nanoparticles. After 
the excitation of a surface plasmon, its energy is trans- 
ferred on the femtosecond timescale to the electronic en- 
vironment, resulting in the heating of the latter. On a 
much longer timescale of typically a picosecond, the equi- 
libration of the electrons with the phonon heat bath re- 
sults in the decrease of the temperature of the electronic 
system with time. Such a process is not included in our 
theory. 

An experimentally accessible quantity^S is the dif- 
ferential transmission AT/T = (T on —T q)/T q, i.e., the 
normalized difference of transmissions with and without 
the pump laser field. It is related to the absorption cross 
section by means of 

AT 3 

(u,T) = -—^ [a(u,T) - a(c,T amb )] , (69) 

Tamb being the ambient temperature. The relation l(69() 
between transmission and absorption holds provided that 
the reflectivity of the sample can be neglected, which is 
the case in most of the experimental setups. The differ- 
ential transmission can be viewed in two different ways. 
For a fixed time delay between the pump and the probe 
pulses, it is sensitive to the energy provided by the pump 
laser which is transferred to the electronic environment 
via the surface plasmon, and thus to the temperature of 
the heat bath. Alternatively, for a given pump inten- 
sity, increasing the time delay between the pump and 
the probe scans the relaxation process of the electronic 
system as the bath temperature decreases. 

In Fig. [SI we present the differential transmission of 
((69(1 for a sodium nanoparticle of radius a — 2 nm, which 
has a strong surface plasmon resonance around 3.5 eV. 
We see that as the temperature of the electronic system 
increases, AT/T becomes more and more pronounced 
since both, the linewidth and the redshift of the reso- 
nance frequency in ((68(1 increase with temperature. Sim- 
ilar results are observed in experiments, 32,33 where the 
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FIG. 8: Differential transmission AT/T as a function of the 
probe energy hu> for increasing temperatures, resulting from 
increasing pump intensities and a fixed delay between pump 
and probe (or by decreasing the time delay at a fixed pump in- 
tensity). The presented results are for a sodium nanoparticle 
with a radius a — 2 nm. 



amplitude of AT/T is observed to decrease as a function 
of the time delay between the pump and the probe laser 
fields. This is accompanied by the blueshift of the cross- 
ing of the differential transmission curves with the zero 
line as the time delay increases. Moreover, the asymme- 
try of AT/T is observed in the experiments and obtained 
in our calculations. 

Our results could provide a possibility to fit the ex- 
perimental results on metallic nanoparticles excited by a 
pump laser field in order to extract the temperature, and 
thus could assist in analyzing the relaxation process. 

VII. CONCLUSION 



potential. The ratio of spill-out electrons over the total 
number as well as the resulting redshift of the surface 
plasmon frequency scale inversely with the size of the 
nanoparticle. The spill-out-induced redshift was shown, 
by means of a Sommcrfcld expansion, to increase with 
the temperature. 

We have demonstrated that the coupling between the 
electronic center of mass and the relative coordinates re- 
sults in an additional redshift of the surface plasmon fre- 
quency. This effect is of the same order as the redshift 
induced by the spill-out effect and presents a similar size 
and temperature dependence. Thus it has to be taken 
into account in the description of numerical and exper- 
imental results. Our semiclassical theory predicts that 
for the smallest sizes of nanoparticles, the environment- 
induced redshift should exhibit a nonmonotonic behavior 
as a function of the size, as confirmed by numerical cal- 
culations. This is not the case for the redshift caused 
by the spill-out effect, and thus permits to distinguish 
between the two effects. 

Our theory of the thermal broadening of the surface 
plasmon resonance, together with the temperature de- 
pendence of the resonance frequency, qualitatively ex- 
plains the observed differential transmission that one 
measures in time-resolved pump-probe experiments. Our 
findings could open the possibility to analyze relaxation 
processes in excited nanoparticles. 
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The role of the electronic environment on the surface 
plasmon resonance in metallic nanoparticles has been an- 
alyzed in this work. By means of a separation into col- 
lective and relative coordinates for the electronic system, 
we have shown that the coupling to particle-hole excita- 
tions leads to a finite lifetime as well as a frequency shift 
of the collective surface plasmon excitation. The size and 
temperature dependence of the lincwidth of the surface 
plasmon has been investigated by means of a semiclassi- 
cal evaluation within the mean-field approximation, to- 
gether with a low-temperature expansion. In addition 
to the well-known size dependence of the linewidth, we 
have demonstrated that an increase in temperature leads 
to an increasing width of the resonance. The effect of 
finite temperature has been found to be weak, in quali- 
tative agreement with the experimental results. 

We have analyzed the spill-out effect arising from the 
electron density outside the nanoparticle. Our semiclas- 
sical analysis has led to a good agreement with LDA cal- 
culations. In order to achieve this, it was necessary to 
introduce an effective radius of the nanoparticle which 
accounts for the details of the self-consistent mean-field 



APPENDIX A: LOW-TEMPERATURE 
EXPANSION FOR INTEGRALS INVOLVING 
TWO FERMI FUNCTIONS 

In this appendix, we calculate the function T defined 
in (|32p. which involves the integral of two Fermi distri- 
butions. Integrating (|32fl by parts yields 



with 



H{x) 



dx'H(x'), 



(Al) 



(A2) 



the function H being defined in In (SI) , 

we have 

introduced the notation F(e) = [1 — /(e)] f(s — tiu). In 
the low-temperature limit Hlo fceT, we get 



a /3( £ - M ) 



^@(e — frui — [i) 



dF 

-— j+/3 (A3) 

de [©8 + i] 2 r e ,8( s -fiw-j0 4- l] 2 
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which corresponds to two peaks of opposite sign centered 
at e = fi and at e = fi + Hlj. It is therefore helpful 
to expand in I|A1(I the function Ji around e = p and 
e = fi+fko. For low temperatures and for \fku — /i\ > k^T, 
we obtain 

- « (l + #;)-«(#;) (A4) 



7T 



H' 1 



where if' denotes the derivative of if. This treatment 
is appropriate for all values of u>, except in a range of 
order k^T around /i. There, the nonanalyticity of TL 
has to be properly accounted foi— since for x near 1, 
H{x) ~ (16/15)(x - 1) 5/2 . At So; = M> an additional 
term thus appears, and 



^( M)M )~W(2) 



6 V ftw 



if '(2) - C 



So; 



with 



C 



da;- 



o (e* + iy 



3.07. 



5/2 



(A5) 



(A6) 



In order to pursue the evaluation of (|A4|) . we need 
to determine the chemical potential \i. Since the DOS 
of QH()[I. once summed over I and m yields^ the three- 
dimensional bulk DOS proportional to y/e, we can use 
the standard Sommerfeld expression for the chemical po- 
tential of free fermions4£ We thus get for \tvjj — fi\ > k^T 



Fiji, hw) 



£F 



1 

2 



if' 1 



(A7) 



where Tp is the Fermi temperature. Inserting the result 
(|A7|) into i|31|) . and using the expression (|33[1 for if, we 
finally obtain J21J. 
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